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Abstract 

This paper describes real-time attitude determination results for the Solar, Anomalous, and Magnetospheric Particle 
Explorer (SAMPEX), a gyroless spacecraft, using a Kalman filter/Euler equation approach denoted the Real-Time 
Sequential Filter (RTSF). The RTSF is an extended Kalman filter whose state vector includes the attitude quaternion and 
corrections to the rates, which are modeled as Markov processes with small time constants. The rate corrections impart 
a significant robustness to the RTSF against errors in modeling the environmental and control torques, as well as errors 
in the initial attitude and rates, while maintaining a small state vector. 

SAMPEX flight data from various mission phases are used to demonstrate the robustness of the RTSF against a priori 
attitude and rate errors of up to 90 deg and 0.5 deg/sec, respectively, as well as a sensitivity of 0.0003 deg/sec in 
estimating rate corrections in torque computations. In contrast, it is shown that the RTSF attitude estimates without the rate 
corrections can degrade rapidly. RTSF advantages over single-frame attitude determination algorithms are also 
demonstrated through (1) substantial improvements in attitude solutions during Sun-magnetic field coalignment and 
(2) magnetic-field-only attitude and rate estimation during the spacecraft's Sun-acquisition mode. 

A robust magnetometer-only attitude-and-rate determination method is also developed to provide for the contingency when 
both Sun data as well as a priori knowledge of the spacecraft state are unavailable. This method includes a deterministic 
algorithm used to initialize the RTSF with coarse estimates of the spacecraft attitude and rates. The combined algorithm 
has been found effective, yielding accuracies of 1.5 deg in attitude and 0.01 deg/sec in the rates and convergence times 
as little as 400 sec. 

1. Introduction 

A crucial aspect of an attitude Kalman filter is that the attitude quaternion q and the covariance matrix must be 
propagated between measurements using the spacecraft’s angular velocity (also referred to here as the rates), «. Thus, 
q is propagated via 


where 


dq _ 1 
dt ~ 2 
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and o x , o> y , and o z are the components of <■> along the spacecraft’s body axes. Note that matrices, including vectors 
(column matrices), are denoted in boldface in this paper, matrix transposes are denoted by the superscript T, and 
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Equation (3) defines a cross product- to-matrix equivalence that will be used elsewhere in this paper. Attitude here 
denotes the orientation of the spacecraft body frame with respect to the geocentric inertial frame (GCI), a describes 
the rate of change of the attitude, and all vectors are resolved along the body axes unless specified otherwise. 


For a gyro-based spacecraft, gyros accurately measure o, and accurate results can be obtained by integrating 
Equation (1) alone. For a gyroless spacecraft such as the Solar, Anomalous, and Magnetospheric Particle Explorer 
(SAMPEX), however, we must propagate the rates via the following Euler equation (References 1 and 2): 


— = N - w x L 
dt 


( 4 ) 


where Nfgf is the sum of the external torques acting on the spacecraft, and L is the total angular momentum of die 
spacecraft* If I is the inertia tensor of the rigid part of the spacecraft and h denotes angular momentum contribution 
from rotating parts such as momentum wheels, L is given by 


Equations (1) and (4) generally pose a difficult initial-value problem, because significant errors in u can arise from 
spacecraft and torque models, as well as the initial values of q and <■> used in the integration. 


Of course, we may avoid rate propagation altogether by using only single-frame algorithms such as TRIAD (algebraic 
method in Reference 1) and QUEST (Reference 3), which determine the attitude using at least two simultaneous 
measurements. But these methods do not provide accurate solutions when all of the observed vectors are nearly collinear. 
A Kalman filter is superior in such a situation, as will be demonstrated later, since it can use propagation to estimate 
the unobservable elements of the state. Additionally, modeling the spacecraft’s dynamics is advantageous because it 
accounts for other physical phenomena such as nutational frequencies. 


A Kalman filter scheme that corrects for the errors in o> is thus desirable. This paper presents results from such a filter 
developed for the PC-based SAMPEX real-time attitude determination system (Reference 4), denoted the Real-Time 
Sequential Filter (RTSF). Only pertinent aspects of the SAMPEX RTSF are discussed here; complete details, including 
system aspects, mathematics, and models for the environmental torques, are given in References 5—7. 

The basic features of the SAMPEX RTSF are as follows. Let « denote the rates generated through Equation (4), 
and let denote the true rates. The errors are the difference of the two, and these errors are modeled as being of 
two types: a zero-mean white noise vector, q tt , and a systematic error, i.e., 

u = u + b , + n ® 

prop oiu out ■« 

In addition to the quaternion, we wish to estimate the rate errors b using the following Kalman state vector: 

X = [q T b (7) 

together with the following general principles: 

• Model N m as accurately as possible. 

• Model noise terms such as by treating them as a noisy background arising from torque uncertainties and 
accounting for them statistically in the Kalman filter formulation. 

• Estimate b optimally by using the sensor residuals, and propagate b between measurements using a suitable 

dynamics model. 


To formulate a dynamics model for b , we first note the following aspect of Equation (4): It is not reasonable to treat 
b as constant because even constant errors in N UI do not necessarily generate constant b In fact, in view of the 
sporadic nature of the control torques as well as the attitude dependence of the environmental torques, it is uncertain 
whether a simple dynamic model can be developed for b. We avoid this issue by arguing that we are interested not in 
the origin of b , but in its instantaneous value that would reconcile the differences between the propagated and observed 
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values of q. Consider now propagation between measurement times t k and t k ^ and use the following notation: filter 
estimates are denoted by carets; estimates before and after updates (via measurements) are denoted by minus signs and 
plus signs, respectively; and the times are denoted by subscripts k and *+ 1 . In this notation, d)-) = o> (n and 
we seek to correct the rates using * * ' 

**( + ) = &*(-)- * t ( + ) ( 8 ) 

Propagation of q, « and the covariance matrix, P, to t is then done using 6 t ( + ). Since Equation (8) corrects 
tne rates at^ t k , it follows that must reflect only errors that accumulated during the intervening period, 

Kalman filter updates, however, are of the form: 

^i( + ) = + w 

where Abj. denotes the correction estimated by the filter. It follows that ^,( + ) (and Ab’ kfl ) will denote the errors 

m eu j* ^ = ° '^ 1JS * s ac hi ev ed in the RTSF by specifying that 6 decays exponentially with a time constant, 
x, of the order of At. That is, we model the dynamics of b through the first-order Markov model: 

db _ lt 

* ■ b + ( 10 ) 

where i| 6 is another zero-mean white noise vector. 

"Hie above Markov model also assists in dining the filter through the following useful feature: The diagonal elements 
of P _ corresponding to b would then converge to a constant value that depends on x and the statistical properties of 
Forsimphcity, take b and i\ b as scalars, and define Q, through <t) 6 (0 Tj 6 (r') 7 > = where is 

the Dirac delta function. Then p, the Kalman filter’s covariance matrix element for b, converges to p_ given by 


Let us now suppose that the < |tj u | 2 > is known by examining the model uncertainties in Equation (4). We argue that 
P~ ,s so 0 “ us or ^ er °f magnitude since b cannot be known to a greater degree of precision than w itself. Thus, we 
use p_ and x together to fix the numerical value of Q in Equation (11), and then use Q in the Kalman filter 
computations. 

This formulation for estimating b has many similarities to gyro-bias estimation (References 8-10), where b denotes the 
gyro drift-rate biases and «(-) donates uncorrected gyro rates. In view of the slow variation of these biases over time 
(see, e.g.. Reference 1 1), it is then possible to approximate their dynamics by Equation (10) using a large value for x 
There is one important difference, though, between the two situations. Gyros accurately measure w, and the leading 
errors in the gyro rates are indeed removed by estimating essentially constant biases. Thus, for example, 6 (*) is 
subtracted from both «*(-) and and before their mean is used to propagate q via Equation (1). Consequently, 

e increments, An, to n provided by the Kalman updates take 6 progressively closer to that constant. 

Th f r ® St ° f this P aper 1S organized ^ follows. Sections 2 and 3 summarize details of the spacecraft, the data used here 
and the theory. Sections 4 i and 5 are devoted to tuning, accuracy, and robustness of the filter. Section 6 demonstrates 
the advantages of the RTSF over single-frame solutions by (1) showing its stability when Sun and magnetic field vectors 

T Pai ? ! T i (2) delenmrUng attitude and rates usin S only magnetic field data. Section 7 shows that combining 
die RTSF with the deterministic algorithm of References 12 and 13 yields a robust magnetometer-only attitude-and-rate 
determination scheme. Section 8 summarizes the conclusions. 

2. Description of SAMPEX and Data Characteristics 

^ PEX °f the 5ma11 Ex P Iorer satel]ltes and is designed to study elemental and isotopic composition of 

au C °/ miC ° rigin - 11 has a 550x675 ^ orbit with “ 82-deg inclination. SAMPEX nominally 

is Sun-pointing and has a rate of 1 rotation per orbit (RPO) about the spacecraft-to-Sun vector. The attitude accuracy 
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requirement of 2 deg is achieved using a fine Sun sensor (FSS) and a three-axis magnetometer (TAM). Tht control 
hardware consists of a momentum wheel and a magnetic torquer assembly (MTA). 


The nominal body y-axis is the Sun vector, and the nominal x- and z -body axes are defined by the Sun/orbital frame 

of Reference 14 For our purposes, it suffices to note the following. 

. The body y-axis is the pitch axis and also the FSS boresight, and the angular delations of this axis from the 

Sun vector are defined as the roll and yaw angles, which are measured by die FSS. 

• The pitch angle denotes angular errors about the y-axis and is measured by the TAM. • 

. Pitch P control is derived through the wheel whose axis coincides with the body y-axis, whereas the MTA 

primarily provides roll and yaw control, and angular momentum dumping. 

. During the nominal 1-RPO mode, roll, pitch, and yaw angles are all 0, and « = (0, 0.06,0) deg/sec. 


SAMPEX telemetry datasets of Table 1 are used here and, with the exception 
SIM722 was constructed (Reference 7) using SAMPEX ephemerides and the 
calibrating the accuracy of the b estimates. 


of SIM722, they all contain inflight data, 
wheel data from PB722 and is useful in 


Table 1. Highlights of the SAMPEX Telemetry Datasets Used for Evaluating the RTSF 


Mann® of Dataset 

Features 

PB705 

Transition from Sun-acquisition mode to nominal 1-RPO mode on the day of launch, 7/4/92 

PASS53 

Nominal 1-RPO mode data of 7/8/92 

SIM722 

Simulated data for duration of eclipse (Earth shadow) on 7/21/92 

PB825 

Near-coalignment of Sun and magnetic field vectors on 8/23/92 

PB722 

Inflight data for the eclipse of SIM722 


No significant differences were noticed between QUEST and TRIAD attitude solutions. The TRIAD attitude solutions 
were differenced to provide rates, and the single-frame results are often taken as the truth models accurate to. 0.5 d g 
for roll/yaw, 1.5 deg for pitch, and 0.5 deg/sec for rates. 


3. Theoretical Aspects of the RTSF 


Definition of Errors: 

The Kalman filter formalism here follows the scheme of Reference 15; i.e., state errors before and after an update are 
considered and the differential equation for P is derived by demanding optimality. The following notation is used m 
addition to that in Section 1 A generally denote an orthogonal matrix; A(q) is then the orthogonal matrix parameterized 
^ S^ntion multiplication, denoted by ® . is defined here in the reverse order of the corresponding 

attitude matrices; i.e., q ■ q l ®q 2 is equivalent to A(f) = A{q}A{qJ. 

The attitude error is linearized, and is taken as a vector of three small independent Euler angles, a ■ (a,, a y a}, needed 
to rotate the true body frame onto the estimated body frame. That is, if 


6q 


•I 1 


( 12 ) 


then A(4) = A(dq)A(q nt ) ■ The error x in the state vector X of Equation (7) is then 

x ■ [o r Ah r ] r 


(13) 
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The 6x6 covariance matrix P is then defined as P ■ <xx T >. If X 
following relationships hold: 


I T T 17 * 

Qnu tih* tnie state vector, the 


4(±) = 4tnu® *9(±) 


^(*) = ~ AA(±) 

The sign convention for AA ensures that, from Equations (6) and (8), the rate estimates are of the following form: 


(14) 


6 - <■>„* * Ab + ij 


(15) 


Let 


** = [a* 7- AA* r ] r 

be the errors estimated by the filter. X is then updated through 

4 (+) = 4(-)®(6qr l 

k+) ■ H-) + AA* 

The relation between the state errors before and after an update follows as *(+) = *(-) - x' . 


(16) 

(17) 


Propagation From t k to t k ^ : 


$ is propagated readily after ignoring in Equation (9), i.e., 

^*i(~) = ft( + ) «xp 



Regarding the rates, & is first updated via Equation (8), and L is updated using Equation (5): £ k (+) =/«(+) + A 
Here A is the wheel momentum given by A = [0 /**, 0] r , where is the moment of inertia of the wheel about 

its axis and is the wheel speed. «*(+) and L k (+) are then used to numerically integrate Equations (1) and (4) to 
obtain $**,(-) and £^,(-), along with the differential equation for P to be given shortly. The uncorrected rates at 
f **i are then obtained by inverting Equation (5), i.e., 6 i+1 (-) = -A j. 


While propagating L, N m is modeled as the sum of four external torques: magnetic control ( N ), gravity gradient, 
aerodynamic, and radiation pressure. (The wheel torque is implicit through h k ^ in the above propagation scheme.) The 
environmental torques are computed in the RTSF using a detailed spacecraft model as well as accurate algorithms; in 
particular, the aerodynamic torque is computed using the FREEMAC algorithm (Reference 16), which generates velocity- 
dependent drag coefficients. We note, however, that the noise in SAMPEX wheel speed obscures the effects of the 
environmental torques. Thus, only N is noticeable; this is given by N = pxfi, where p is the dipole moment 
of the MTA and B is the magnetic field. 


The differential equation for P can be derived by noting that the error in £ after propagation comes from two sources: 
(1) the error before propagation and (2) the random vectors, tj o and ij 4 , which are ignored in propagation. The 
linearized propagation equations for the errors follow as 


dx 

*[«x] 

Ax3 



dt " 

°3*3 


X + 

.V 


» F(a)x + q 


(18) 


Equation (18) defines F and i). Introducing Q, the diagonal spectral density matrix of tj, through 
<T l(0 'nftj > = Q bit- 1 1 ) it can then be shown (Reference 15) that the covariance matrix evolves according to 

? = FP * PF T * Q (19) 

at 

Equation (19) is integrated numerically in the RTSF to obtain P k , x {-). 


485 


Updates Using Sensor Measurements : 

Let V and V be the measured (body frame) and reference (GO) vectors corresponding to a particular sensor and let 
V V J C-)) C tL residual y is defined as y = V* - V Linearizing the attitude errors about yields 

5 y - *,(-) ♦ AV, (2 °> 

where H ■ f-[V xl 0, ,1 and A V B is random measurement error. The foUowing optimahty requirements can then be 
obti K uTkalman gafn and covariance-update equations: (1) if the error is unbiased before ameasme^t 
-)> = 0), it should remain so after the update (<*(+)> = 0), and (2) the statical average of the total error, 

<x(+) r x(+)>, should be a minimum: 

K = P(-)H T [HPi-)H T ♦ *]-‘ (21) 

*♦) = [i 6x6 -**]*-> 

where R ■ <AV_ AV» r >. The correction x* is then given by 

( 22 ) 

x* = Ky 

for use in Equations (16) and (17). The covariance update in Equation (21) was chosen over more numerically: stable 
algorithms due to memory and speed considerations on personal computers. To partially mitigate the risk s, *+) 
symmetrized after the update. No adverse effects have been noticed to date. 

For SAMPEX no special treatment is necessary to construct R for the TAM, since independent magnetic field 
measurements are madTalong each body axis. The FSS, however, measures only twoind^endent angles, ando y 
corresponding components of y are used in the computations. Details are given in Reference 5. 


4. Tuning and Accuracy of the Rate Corrections 


The telemetry digitization of the FSS is 0.5 deg; this value was used to construct Jt for FSS me^uremen ± 
TAM uncertainties in the reference magnetic field calculations are larger than the digitization errors of 0.3 mG, 
the residuals were examined, a TAM noise of 3 mG was used to construct R 

A distinctive feature of telemetered SAMPEX data is the large wheel noise (-1 rad/sec). This no ^ £ 

. j in* 2 N m far in excess of the maximum values of the environmental torques ( 10 N m), 
r^of ata. L7xio“ IXcTL . The pitch covuhauce f n thus gmws by about 3xU>*(W rad’ between 
n ■ the rate 0 f chance of P due to the process noise (Equations (19)), Q u was correspondmgly 

- ta t ate substantial (as during Sun 

acquisition), the same value was also assigned for <2 n and Q n . 

The discussion of Equations (8-11) then fixes the bottom three elements of <?, which correspond to i\ b . Thus, p m m 
Equation (11) was chosen as 3x10'* rad’/sec’, so that 

, «*12l 2*! , | . 4. 5, 6 

“ r sec 3 

T is chosen equal to the telemetry period: 0.5 sec for real-time data and 5.0 sec for playback data. As can be seen from 
Figure 1, the rate-error elements of P converge to the above value of p m very well. 

Retire 2 clarifies the discussion of Equations (9) and (10), by examining the dependence of t on the up<kte penwl 
fofl mtTs Tis updated via EquVtton (17) only after every sec, as opposed to the telemetry penod of 5 sec here. 
(S is alwavs uodated^every 5 sec.) Here "Sun angle" is the angle between the predicted and measured Sun vectors, thus, 
ftr a tJsS meL^meX larg^ Sun angles indicate larger propagation errors. Plots a and b show that when 
x = r t 5, the mean Sun angle is about 0.25 deg and rate corrections ofthe order of 10 deg^o^ are wttmated. In 

vpd 
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ause 6 does not decay to zero during the time between the updates, so that the updates, £(+), in Equation (9) also 
include past rate errors. The divergences are eliminated by choosing f = 25 sec, as shown in plot fHowever since 
6 is now updated infrequently, the Sun angles in plot f are generally T^ger than those in plot l 

6 CStimates V 35 ascertained b y usin 8 simulated data. These results are shown in Figure 3. In plot a 

dl wlTn^u to l t0 fluCtUati0n ^ in infli S ht wheel data ^ed to generate SIM722. The simulated 
data were input to the RTSF but mth a constant wheel speed equal to the value at 0 sec , and we see that the RTSF’s co 

c f^ a ^ S . generaIIy refle, : t true “ ean V£dues - However, significant differences occur around 100 and 450 sec, where 
g es tn cause changes in the true w This information is not available to the RTSF and leads to significant 
Sun angies, together with nonzero values for B y around these times in plot b. We see from plots a and b that the signs 
ese B y estimates agree with Equation (6). Quantitatively also, the B estimates are reasonable; for example, the 
corrections total about 0.01 deg/sec during 0-200 sec, which compares well with the rate differences in plot a. 

nL^ouT ^ tUmng rCSUltS C ° nfmn * e argUmentS of Section 1 md «te corrections B are accurate to 

5. Robustness of the RTSF Against A Priori Errors and Torque Errors 

^Ift F WaS f eV n 1Uated CXtensiv J f y with PB705 data where SAMPEX is in Sun acquisition soon after launch. The 

atouTofdeVZ Lf ^ raP y_W > u C ° nStant at about 0 3 de 8 /sec and and sinusoidal with amplitudes of 
S ' Tn 2^ o fransitiomng to the nominal 1-RPO mission mode. (The transition is clearly visible in Figure 

■ us, e un-acquisition part of PB705 is a rapidly varying situation and provides a stringent test of the RTSF’s 

uSSSd , "ot LTiT^m RTSF SUned Wi,h What ““ body fra™ 


Figure 4 shows the convergence of the RTSF results to single-frame solutions using both FSS and TAM data Although 
the a pnon errors are large— (-27, -91,21) deg in attitude and (-0.5, 0.3, 0.4) deg/sec in the rates— the RTSF’s estimates 

convergence 250111 3 °° ^ ^ ° f ** 3 Pri ° ri erT ° rS Clearly show up in p,ot c ’ where ^ is extrem ely large during 

Figures 5 and 6 i show another useful aspect of the RTSF: the ability to compensate for torque errors Figure 5 shows 
results with PB705 during the transition to 1-RPO mode, with the TRIAD rates serving a? the truth m<Ll. Hereo 

input to Te RTSF^ ^ ^ ^ ° f 006 deg/SCC P3rtly due t0 3 substantial wheel torque. The wheel speeds 

SI t Ti : ever ’^, ere dellberatel y ven wrong sign during this run. Dius, around 200 sec, the RTSF’s « 

coS ™ \ mCre , aSeS - ^ glves rise t0 Iar « e Sun “ d significant B y are estimated (plot b) that eventually 

y 'u ag r needed t0 C0ITeCt “ iS fmally eHminated in the 1-RPO mode where « is nearly constant 

SSSSXTdS? dUri " 8 ““ ° f S' 0 " 1 "' '° dM “' conversion factor 

l Pre f ntS 1 c UltS With PASS53 ^ where rate e mors of about 1 deg/hour appear in phase with MTA activity 
Shtauld m ^ 5 ’ d3ta hCre 316 unmodified > and the RTSF results suggest that the MTA needs to be 

SSL 7 PrCSentS reSU,tS witb PASS53 data highlighting the usefulness of estimating b even in a slowly varying 
“ H f. re P ro pagation of u via Equation (4) was completely omitted, and the RTSF was run using zero initial 

data where \h ° >y eStUn3te quickly converges to the 1-RPO value. In contrast. Figure 8 shows the results with PB705 
date where the reverse situation holds: Ifre a priori errors were small, « was propagated, but b was not estimated 
ge rate errors (which we ascribe to errors) accumulate over time and eventually degrade the attitude accuracy. 

6. RTSF Advantages Over Single-Frame Algorithms 

£e ro^^ ti0n ’ a “Tf" fil f yieldS estimates Via propagation of ^ for stat o vector elements that 
" K b l K m ! demonstrated here for two situations: (1) near-coalignment of the Sun and magnetic field 

dS ^riTh^F^ 6 PU ?- “f ® 1S ne3r,y unobservable ( 2 ) attitude and rate estimation using only magnetic field 
ta. (Note that Figure 3 is also a magnetic-field-only situation, since SIM722 models an eclipse 8 ) 
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vectors whereby the pitch angle is nearly unobservable and (2) attitude and rate estimation using only magnetic field 
data. (Note that Figure 3 is also a magnetic-field-only situation, since SIM722 models an eclipse.) 

Fitrure 9 presents the former case, with the top plot showing peaks of nearly 25 deg in both QUEST and TRIAD 
™Ens for the pitch angle when the Sun-magnetic field angle falls below 5 deg. This occurs despite nominal values 
£5 control data (not shown here) during this period. In contrast, the RTSF pitch estimates vary very httle 
over the timespan with these estimates now being generated using propagation alone. The absence of pitch observabi ity 
~SrTS pitch covariance P n , whose growth due to <? in Equation (19) is not offset by measurements, 
first rises Steeply in plot b amfthen falls when the TAM data again prov.de useful information. 

Fitrure 10 presents RTSF results using PB705 data but without using the available FSS date. The starting conditions 
here were ET Sal conditions Jin Figure 4, but now attitude information is obtained through oriy he , -gnetm 
field B Although a single measurement of B can yield only information regarding rotations perpendicu ar o t 
changes directiof sufficiently during an orbit, so that all three attitude angles are observable want mffic endlong 
span of data. This is indeed seen in Figure 10, where the RTSF’s yaw angle estimate converge after about 1000 sec^ 
After convergence, the RTSF attitude and rate estimates remain within 1.5 deg and 0.01 deg/sec, respec y, 
single-frame estimates (not apparent here). We attribute this difference to magnetic field model uncertainties smee RTSF 
attitude accuracies of 0.1 deg were obtained using simulated data (Reference 7). 

7. Robust TAM-Only Attitude and Rate Determination 

Motivated by the successful T AM-only attitude-and-rate estimation by the RTSF 3 ®“ d 10) ’ we ha ^® 

a robust TAM-only method that overcomes the following shortcomings of the RTSF: (1) Convergence 
LfirmSs ifVot guaranteed since it is an extended Kalman filter, and (2) the TAM-only convergence is slow about 
STJS hr Figure^” 1 These difficulties are overcome here by initializing the RTSF with the solution from the TAM- 
only deterministic algorithm (DA) of References 12 and 13. 

The DA considers the following two time derivatives of the magnetic field the one in the J 

the one in the body frame, B A . (The notation used here conforms with that of References 12 and 13.) They are related 
formally by the vector equation: 

(23) 

B x = B a + <**B 

Noting that |B| is invariant under orthogonal transformations, two special orthogonal axes peipendicular to 
B <£efe“ e 12) can be chosen such that to projection, o, and te, of «e along these «es are restarted to he on 
a circle. The DA then extracts the radius and center of this circle from TAM data. 

The constraint on <o 2 and o> 3 can be parameterized by an angle * so that « is described by two unknown parameters: 
the projection of <o onto B, and 4>. These are then determined using « and die second time de " vatives of 
as evaluated in the two frames. From Equation (4), d> is expressed as a function of w and N ut , yieldmg 

(24) 


+ + = 0 

. the vect ors A (4>) it = 0, 1, 2 , are defined in Reference 13. Projecting Equation (24) onto the plane 
perpendicular to B yields two transcendental equations for * and o>„ generally yielding at least two solutions_ The 
Lbiguity in the solutions can be removed if N M can be neglected, so that we can demand that the reference frame 
components of the angular momentum L be constant. For SAMPEX, in particular, this requirement holds during 

eclipses when there is no magnetic control. 

The efficacy of the DA calculations was tested using eclipse data from PB722, and evaluating B R and Bousing 100 sec 
Figure 1 1 preset attitude Shd ode results fro,, this emulation. tmtng a 2-3-2 Euler seque.ee 
iT^terSTthe attitude. The advantage of this pammeterization is Out. the third Euler angle three, ly reflects the 
1-RPO rate of the spacecraft, whereas the other two angles are very nearly constant for small nutattonal ampl tudes. 
We see that although up to four solutions appear toward the middle of the timespan, only two solutions appear owar 
r e'reu^ SLtding dm consuncy ofi in the refetence frame dte. unmnblguously selecm dte coaect solutton 
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Why We USUally have ° nly tW0 solutions - k is useful to approximate L ~ h, which holds very well for 


«(*,<*!) -q2(*) + Ox Q? 


(25) 


where QoO®) ■ x *]> Qj |£| ■ [BxA], and co i denotes the projection of the o onto the plane 

perpendicular to B. This linearizes Equation (24) and yields: 


AoW + OxAjW-O (26) 

One can then exclude o, by projecting these equations onto the vector C’ 3 \B\ - BxA,(4>) . This leads to a quadratic 
equation with respect to x = tan(4>/2) , analogous to that in Reference 12 for the constant- o limit. (The quadratic 
equation turns into Equation (3-11) in Reference 12 if one sets h = 0 .) 


The DA yields multiplesolutions when B is approximately antiparallel to the roll axis. Since h for SAMPEX is directed 
along the pitch axis, is approximately antiparallel to B for any value of 4> . Contributions to o2(®) cannot then 

be neglected m Bxia, and this gives rise to the multiple solutions, at about 600 sec in Figure 11. 

Figure 12 compares the RTSF roll angle results obtained by initializing the filter with two different schemes: (1) zero 
nutial conditions and (2) with the correct DA solution from Figure 1 1. It is evident that using the DA for initialization 
significantly improves the attitude accuracy by reducing the convergence time of the RTSF. Note that results with both 
initial conditions reflect oscillations with the spacecraft’s nutational period (=120 sec). To illustrate the importance of 
the initial conditions, we note that the RTSF solutions diverged when the spurious solution of Figure 11 was used to 
initialize the filter. In this context the successful convergence with zero initial conditions is also noteworthy. (The large 
nutational amplitudes for zero initial conditions are eventually damped by the RTSF’s rate corrections.) 

8. Conclusions 


We find that the SAMPEX RTSF’s method of estimating the rate errors b that accumulated during the period between 
rate updates yields accurate results while minimizing the computational load. The tuning scheme, which exploited the 
relationship between the converged covariances, noise characteristics of SAMPEX data, and the Markov time constant, 
was simple but effective, and resulted in estimates of b accurate to 0.0003 deg/sec (1 deg/hour). Consequently 
estimating b significantly enhanced the performance of the SAMPEX RTSF by providing robustness against 

• Large a priori attitude and rate errors 

• Errors in the control torques 

In fact, we see that rates were generated accurately in slowly varying situations even when the Euler equation is not 
used. 


The RTSF’s advantages over single-frame algorithms were demonstrated through 

Stability of pitch angle estimates during Sun-magnetic field near-coalignment, where both QUEST and TRIAD 
solutions were incorrect by nearly 25 deg 

• Magnetic-field-only attitude and rate determination to within 1.5 deg and 0.01 deg/sec of the single-frame 
solutions using zero initial conditions 


The deterministic magnetic-field-only algorithm significantly enhanced the robustness and accuracy of the filter by 
generating coarse estimates of the a priori attitude and rates. This results in an important contingency algorithm for 
spacecraft like SAMPEX which have no sensor redundancy. 
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SUN ANGLE (DKG) 



Figures 2a-2d. RTSF Results Using PB705 Data Showing That Tuning is Optimal 
When Both x and t upd Equal the Telemetry Period of 5 sec 
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Figures 2e and 2f. RTSF Results Using PB705 Data Showing That Tuning is Optimal 
When Both x and t upa Equal the Telemetry Period of 5 Sec 
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Figures 3a and 3b. Results Using SIM722 Data Showing Accuracy of the Rate Corrections 
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Figures 4a and 4b. Results Using PB705 Data Showing Convergence of RTSF Estimates 
to Single-Frame Results Despite Large a Priori Errors 


















Figure 7. Results Using PASS53 Data 
Showing Robustness of RTSF When Rate 
Propagation Is Suppressed 


Figure 8. Results Using PB705 Data 
Showing Significant RTSF Errors When Rate 
Corrections Are Not Estimated 
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Figure 9a and 9b. Results Using PB825 Data Showing Robustness of RTSF During 
Near-Coalignment of Sun and Magnetic Field Vectors 




Figures 10a and 10b. Convergence of RTSF Estimates for PB705 Data When Only 
TAM Data and Large a Priori Errors Are Used 
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Figures 11a and 11b. T AM-Only Attitude and Rate Solutions Generated by the Deterministic 

Method During an Eclipse in PB722 Data 



Figure 11c. TAM-Only Attitude and Rate 
Solutions Generated by the Deterministic 
Method During an Eclipse in PB722 Data 



Figure 12. RTSF TAM-Only Results for 
PB722 Data Showing Faster Convergence 
Using the Correct Solutions of Figure 1 1 
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